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Abstract 

We calculate the full asymptotic expansion of boundary blow-up so- 
lutions (see equation ([T} below), for any nonlinearity /. Our approach 
enables us to state sharp qualitative results regarding uniqueness and ra- 
dial symmetry of solutions, as well as a characterization of nonlinearities 
for which the blow-up rate is universal. Lastly, we study in more detail 
the standard nonlinearities f{u) = uF , p > 1. 

1 Introduction 

Let B denote the unit baU of R^, iV > 1 and let / e C^(R). We study 
the equation 



(1) 



Aw = /(u) in B, 

u — +00 on dB, 



where the boundary condition is understood in the sense that 
lim u{x) — +00 for all xo G dB 

X — ^XQ B 

and where / is assumed to be positive at infinity, in the sense that 

(2) 3 a e R s.t. /(a) > and f{t) > for t > a. 

A function u satisfying ([TJ is called a boundary blow-up solution or simply 
a large solution. Existence of a solution of IT} is equivalent to the so-called 
Keller-Osserman condition : 

(3) J < where F{t) = j f{s) ds. 



1 
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For a proof of this fact, see the seminal works of J.B. Keller [8] and R. 
Osserman [11] for the case of monotone /, as well as [6] for the general 
case. From here on, we always assume that Q holds. 

Our goal here is to study asymptotics, uniqueness and symmetry prop- 
erties of solutions. Our approach improves known results in at least two 
directions : firstly, aside from the necessary condition ([3]), we need not 
make any additional assumption on / to obtain the sharp asymptotics 
of solutions. Secondly, we obtain the complete asymptotic expansion of 
solutions to all orders. Here is a summary of our findings. 

Theorem 1.1 Let f £ C^(R) and assume ((2|, ((3| hold. Consider two 
solutions ui , U2 of m . Then, 

lim ui(a::) — U2{x) = 0, for all xq € dB. 

More precisely, there exists a constant C — C{u\,U2, N , F) > 0, such that 
for all X £ B, 

(4) \u^(x) - U2{x)\ < C dt. 



U2(x) 



In addition, 

(5) \F{m)-F{n2)\&L'-{B). 

Estimates on the gradient of solutions can be obtained for a restricted 
class of nonlinearities, namely 

Theorem 1.2 Let f G C^(R) and assume and Q hold. A ssume m 
addition that f is increasing up to a linear perturbation i. e. there exists 
an increasing function f and a constant K such that 

(6) f{t) = f{t) - Kt, for all t G R. 
Consider two solutions u\ and U2 of ([1]). Then, 

(7) \V{u,^U2)\eL^{B). 

The previous theorems can be used to study qualitative properties of 
solutions, such as uniqueness and symmetry. We begin with the question 
of uniqueness of solutions of ([1]). The following conjecture is due to P.J. 
McKenna ([5]). 

Conjecture 1.3 ([S]) Let N > 1, Q, a smoothly bounded domain ofW^ 
and f £ C^(R) a function such that ((2| and ((3| hold. Assume in addition 
that the function f defined by 

(8) f{t) = f{t) - Xit, for allteR 

is increasing, where Ai = Ai(— A;f2) > denotes the principal eigenvalue 
of the Laplace operator with homogeneous Dirichlet boundary condition. 
Then, there exists a unique large solution of 

Au — f{u) in n, 

u — +00 on dQ. 
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As a direct consequence of Theorem 11.11 we prove Conjecture 11.31 in the 
case n = B. 

Corollary 1.4 Let f G C^(R) and assume ((U and ^ hold. A ssume m 
addition that that the function f defined by ((Sj is nondecreasing. Then, 
there exists a unique large solution of ([l}. 

Remark 1.5 Many uniqueness theorems have been established in the lit- 
erature (see e.g. the survey [T]^, and they hold for a general class of 
bounded domains Q. However, in all of these results, additional assump- 
tions on f are needed, such as convexity. 

Proof of Corollarv ll.41 Assume first that / is nondecreasing. Let ui, U2 
denote two large solutions, ft suffices to prove that ui < U2. Assume this 
is not the case and \et u = {x £ B : w{x) > 0} 7^ 0, where w = u\ — U2. 
Working if necessary on a connected component of w, we may always 
assume that a; is connected. Using Theorem ll.il we see that w solves the 
equation 

' a™ = f{ui) - f{u2) > 

: on duj. 

By the Maximum Principle, u; < in tj, a contradiction. 

Assume now that we only have /' > — Ai. Let > denote an 
eigenfunction associated to Ai and let a = w/ipi, where w = u\ ~ U2 
denotes the difference of two solutions. Assume again that uj = {x £ B : 
w{x) > 0} 7^ 0. By a standard calculation, 

V • (v3?V(t) = (f{ui) - f{u2) + \l{ui - U2))VP1 > in LJ. 

We claim that cr = on dco, from which the desired contradiction will 
follow. By Q and the well-known estimate ipi > c(l — it suffices to 
show that 

r+°° -iL. dt 
(9) hm ''"^W^W ^0 

^ ^ x^dB 1 - \x\ 

We shall prove later (see Lemma [2.4|l that there exists a radial boundary 
blow-up solution [/ of ([1} such that U2 > U. Since U is radial, it follows 
from (121) that U'{r) > for r = Ixj close to 1. In particular, 

U" <AU = f{U). 

Multiplying by U' and integrating the above inequality between ro and r 
close to 1, it follows that (f/')'^/2 < F{U) -\- C. Integrating again between 
r and 1, we obtain 

dt 

< 1 - r. 



lu(r) v/2(F + C) 
for r close to 1. So, 

Iu2(x) Wtt) Iu(x) 'Wft) '^^ C 



V^FiUix)) 
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and (|9]) follows. 



□ 



When /' ^ — Ai, uniqueness fails in general. One may ask however 
whether all solutions of ([T]) are radial. H. Brezis made the following 
conjecture. 

Conjecture 1.6 ([4]) Let f £ C^(R) denote a function such that ^ and 
(jSj hold. Then, every solution of ([1} is radially symmetric. 

To our knowledge, the first contribution to the proof of Conjecture 11.61 is 
due to P.J. McKenna, W. Reichel, and W. Walter (see [lO]), using the ad- 
ditional assumption that limt-).+oo /' (t) / ^F{t) = +00. A. Porretta and 
L. Veron then proved the conjecture, assuming only that / is asymptoti- 
cally convex (see [H]). We improve these results as follows. 

Corollary 1.7 Let f G C^(R) and assume that ^ and ((3| hold. Let 
u denote a solution of ([T]). Assume in addition that, up to a linear per- 
turbation, f is increasing (i.e. © holds for some nondecreasing function 
f and some constant K). Then, u is radially symmetric. Furthermore, 



Remark 1.8 In the setting of the classical symmetry result of B. Gidas, 
W. M. Ni and L. Nirenberg (see [7]^, ® is also assumed in order to 
prove symmetry. In the same article, the authors give an example of a 
nonlinearity f failing ([6]) for which there do exist nonradial solutions of the 
equation. In the context of large solutions, we do not have such a counter- 
example. In fact, we expect that none exists i. e. we believe that Conjecture 
holds. But at this stage, we do not even know whether radial symmetry 
continues to hold for simple nonlinearities such as f{u) = u (1 -f sin it). 

Corollary 11.71 is a direct consequence of the moving plane method and 
Theorem O: 

Proof of Corollary 11.71 . Let U denote a radial solution of ([1]). It 
follows from ((2]) that L*' is a nondecreasing function of r = |x| for r close 
to 1~ and ^{r) — >■ -|-cxd as r — > 1^ . By ((7|, we conclude that any solution 
M of lU satisfies ^{x) — >■ +00 as a; — >■ dB, while the tangential part of 
the gradient of it remains bounded. We then apply Theorem 2.1 in [12) . 
□ 

In addition to the relative asymptotic information given by Q, ((5]) 
and ((7|), the exact asymptotic expansion of a solution can be calculated 
to all orders. This is what we explain next. Consider (for simplicity only) 
a radial solution of H]) i.e. a solution of 



with lim^^j- it(r) — -\-oo. We want to think of the second term on the 
left-hand side of (|10|l as a lower order perturbation as r — >■ 1. Multiplying 
the equation by 1; = du/dr and putting the error term on the right-hand 
side, we get 



|i>0^nB\{0}. 



(10) 



d^u N -Idu 



fin), 



dr'^ r dr 



1 d 2 


2 dr 




N -I 



r 



,2 
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Make the change of independent variable u = u{r). Thinking of u as a 
function of the new variable u, we have -r- = ^ tt- = J' :r- and so 

' dr dr du au 



-V 



2 du du r 
In other words, v solves the nonlinear integral equation 



viu) = ^2 (^F{u) j ^dt^ +C =: Af{v), 

where Uo,C are given constants. The above equation turns out to be 
contractive in a suitable Banach space. In particular, it can be solved using 
a standard iterative scheme Vk+i ~ Af{vk)- As we shall demonstrate, each 
Ufc contains (in implicit form) the first k terms in the asymptotic expansion 
of the solution at blow-up. To summarize, we have: 

Theorem 1.9 Let f G C^(R) and assume dgj, ((S} hold. Let Uo G K, 
/ — [Uo, +oo) and let vq be the function defined for u £ L by 



(11) vo{u) = ^/2F{u). 

Consider the Banach space 

X = {v e C(/;R) : 3A/ > such that \v\ < Mvo}, 

endowed with the norm \\v\\ — supi\v/vo\. If the constant Uo is chosen 
sufficiently large, then for some p G (0, 1), there exists a unique solution 
V G B(vo, p) <Z X of the integral equation 



(12) v{u) = ^^2 [f{u)-{N -1) j^^-dt^ 

where r — r{u, v) is given for u £ I , v £ B(vo,p) by 



ue I, 



(13) r(u,v) = 1- - dt 



+ 00 



V 



In addition, v is the limit in X of (vk) defined for k — by (|lip and for 
k>l, by 



and the sequence Vk is asymptotic to v i.e. as u ^ +oo, 

Vk+i{u) = Vk{u) + o{vk{u)) 
and given any k e'N we have 

v{u) = Vk(u) + 0{Vk+l{u) - Vk{u)). 
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Let now u denote any solution of ((TJ and fix ro G (0, 1) such that u{x) > 
Uo for \x\ > ro. For k > 0, define uu for r > ro as the unique solutio'^ 
of 



(15) 



duk , , 

—r- = f fe ) 
dr 



lim Uk{r) — +00, 
where Vk is given by p4|l . Then, 



L 



(r) \Juk(r) VO 



and given any fc G N, uie have 

, , du ( f+°° du\ 

(16) / — = o / — asx-^dB. 

Theorem 11.91 enables one to calculate (implicitly) the asymptotic expan- 
sion of a solution term by term. But how many terms in this expansion 
are singular? This is what we discuss in our last set of results. 

We begin with the simplest class of nonlinearities /, those for which 
only one term in the expansion is singular, namely the function uq defined 
by (|lip and (|15p . It turns out, as A.C. Lazer and P.J. McKenna first 
demonstrated (see [^), that in this case uo{l — d{x)) is the only singular 
term in the asymptotics of any blow-up solution on any smoothly bounded 
domain Q, C and for any dimension N > 1, where d{x) denotes the 
distance of a point a: G fi to the boundary of Q,. In other words, the 
blow-up rate is universal. The question is now to determine for which 
nonlinearities /, this universal blow-up occurs. We characterize these 
nonlinearities as follows: 

Theorem 1.10 Let f2 C R'^ denote a hounded domain satisfying an inner 
and an outer sphere condition at each point of its boundary. Let f G 
C^(R), assume ((2]), ((3| hold and consider the equation 



(17) 

Assume 



Au = f{u) in r2, 

It — +00 on dfl. 



+°° J* V2F ds 



Then, any solution of (|17p satisfies 

(19) \im u{x) - Uo{l ~ d{x)) = 0, 

where d{x) — dist(2;,c'r2) and uo is defined by Hll|) . p5p . 



*Note that JTSj can be solved by quadratures and its solution is unique. Indeed, Vk{u) 
vo(u) = yj2F(u) as ii — >■ +00 and this implies by (jsj that J^°° dt/vk{t) < +oo. 
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(20) 



We also have the following partial converse statement : if 

+°° J* V2F ds 



liminf ^2F{u) 

U — + 



(2f)3/2 



dt > 0, 



then (|19[l always fails. 



Remark 1.11 To our knowledge, p8[l improves upon all known condi- 
tions for H19p to hold (see in particular and [2]j. Despite its unappeal- 
ing technical appearance, (jlSp only uses information on the asymptotics of 
F (in particular, no direct information on f is required). Nonlinearities 
such that F{u) ~ or F{u) ~ u'', p > 4 as u ^ +oo satisfy (fTS)) . For 
F{u) ~ (|20|) /loMs anci so t/ie conclusion H19|) /aiis. 



Remark 1.12 Condition (|18p can 6e weakened to 

ft 



lim v/2-F(wo) 



= 0, 



where uq — Mo(r) is defined by p5|l . Similarly, (|20p can 6e weakened to 



lim inf Y/2F(tii) 



(2F)3/2 



rff > 0. 



As an immediate corollary, we obtain uniqueness on general domains, 
whenever only one singular term appears: 



Corollary 1.13 Assume (|18|l . // in addition, f is nondecreasmg, then 
the solution of (|17|l is unique. 

□ 



Proof. Simply repeat the proof of Corollary 11.41 

More than one term can be present in the asymptotic expansion of u. 
Finding all the (singular) terms in this expansion is of staggering algebraic 
complexity. To illustrate this, we provide the first three terms (in implicit 
form) . 

Proposition 1.14 Let U2 be defined by psp for k — 2. Let also Ri, R2, 
i?3 denote three real-valued functions defined for f/ G R sufficiently large 
by 



( 



du 

TP' 



i?i([/) = (iV- 1) 
R2{U) = (AT - l)x 

iN-l)iy^^^2F 



fV2Fdt 
(2F)3/2 



du. 



2F 



/■+°° ds \ 



dt+ 



V 



+ 



5(A- 1) 



[rv2Fdty 



du 



(2F)3/2 



4 2F 
Then, for all r £ (0, 1), r close to 1, we have 

l-r = Ro{u2{r)) + Ri{u2{r)) + i?2(u2(r))(l + o(l)) 
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For specific nonlinearities, it is possible to invert the above identity. 
This is what we do for f{u) = u^, p > 1: 

Proposition 1.15 Let p > 1 (with 2/(p- 1) ^ and let fiu) = , for 
u > Q. Then, the unique positive solution of ([1} satisfies 



where d{r) = 1 — r for r G (0, 1), and where each Ofc G M depends on N 
and p only. 

Remark 1.16 Proposition \T7T^ was first proved by S. Berhanu and G. 
Porru (see As can be seen from the proof, Provosition M . 1^ remains 
valid for any nonlinearity f such that, for some positive constant c, F(u) = 
cvF'^^ + 0{u) for large values of u (and any solution of the equation). 

Outline of the paper 

1. In the next section, we show that any solution u of ([l]) can be 
squeezed between two radial solutions U and V i.e. the inequal- 
ity [/< it < holds throughout B. 

2. Thanks to this result, we need only find the asymptotics of radial 
solutions to prove Theorem ll.il This is what we do in Section O 

3. To obtain gradient estimates, the squeezing technique is insufficient 
and more work is needed. In Sectional we estimate tangential deriva- 
tives via a standard comparison argument, while we gain control over 
the radial component through a more delicate potential theoretic ar- 
gument. 

4. Section [S] is dedicated to the proof of Theorem ll.9l that is we estab- 
lish an algorithm for computing the asymptotics of solutions to all 
orders. 

5. In Section (6] we characterize nonlinearities for which the blow-up 
rate is universal. 

6. At last. Sections [7] and [8] contain the tedious calculations of the first 
three terms of the asymptotic expansion of u in implicit form for 
general /, and of all terms explicitely for f{u) — u^. 

Notation 

Throughout this paper, the letter C denotes a generic constant, the value 
of which is immaterial. In the last section of the paper, we use the symbol 
Cfc to denote a quantity indexed by an integer k, thought of being "constant 
for fixed fc" , the value of which is again immaterial. 



[2/{p-l)] 




2_ 



as r — >■ 1 , 



fc=0 



2 Ordering solutions 



In this section, we prove that any solution of the equation is bounded 
above and below by radial blow-up solutions. To do so, we impose the 
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following additional condition: g{t) :— f{—t) satisfies ([2]) and 
(21) J -^ = +«,, where G'(t) = /(-t). 



Remark 2.1 Note that (|2ip is not restrictive. Indeed, if u denotes a 
solution of fl} and m = mins u, then u also solves ^ with nonlinearity 
f defined for it G R 6j/ 



/(m) + (m — u) if u < m 
f{u) ifu>m. 



Then, f clearly satisfies (|2ip . 

We now proceed through a series of three lemmas. 

Lemma 2.2 Assume (|2ip holds. For M G E sufficiently large, there 
exists a radial function v£C (B) n C{B) satisfying 

A.V > f{v) m B 

and such that 

V < -M m B. 

Proof. Let g{t) = f{—t) for f G R and let a > Af be a parameter to be 
fixed later on. Since g satisfies ([2]), we may always assume that g{t) > 
for t > M. Let now w denote a solution of 

f-w" = g{w) 
w{0) = a 
w'{0) = 

Claim. There exists an a > M sufficently large such that w{l) > M. 

Note that w is nonincreasing in the set {t : w{t) > M}. We distin- 
guish two cases. 
Case 1. w > M. 

In this case, w is defined on all of R"'". In particular, w{l) > M, as 
desired. 

Case 2. There exists R> such that w{R) = M. 

In this case, since w is nonincreasing in (0, i?), we just need to prove 
that R> 1. To do so, multiply H22p by —w' and integrate between and 
r G (0, R): 

-w' = y'2{G(a)~G{w)), 
where G is an antiderivative of g. Integrate again between and R: 

' -™ dr = R. 



Now, 



R = 



M ^2{G{a) - G{t)) Jo y2(G(^r^GH) 



dt rG-UG(a)/2) 



M V2(G(a) - G{t)) - Jm y/2{Gia) ^ G{t)) " 

G-l(G(a)/2) 



J h 



V^G{t) 
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By (|21|l . we deduce that R> 1 for sufRcently large a. We have just proved 
that w\f^Q > w{l) > M and the claim follows. 

It follows that the function v defined for x G B by v{x) = —w{\x\), is 
the desired subsolution. □ 

Lemma 2.3 Let f £ C(R) and assume ([2| and Q hold. Assume v_ £ C{B) 
satisfies 

Av>f{v) inB. 
Then, there exists a radial large solution V of ^ such that V > v. 

Proof. Let v := A''. Then, v_ and v are respectively a sub and supersolu- 
tion of 



(23) 



Av = f{v) in B, 

v = N on dB, 



provided N is chosen so large that A'' > ||w||l°°(s) and f{N) > 0. Futher- 
more, u < li in B for such values of N. By the method of sub and super- 
solutions (see e.g. Proposition 2.1 in [^), there exists a minimal solution 
Vjv of (|23|l such that A'^ > Vn > v. Note that Vn is radial, as follows from 
the classical symmetry result of Gidas, Ni and Nirenberg (see [7]). Also, 
since Vn is minimal, we have that the sequence (Vn) is nondecreasing 
with respect to A'' (apply e.g. the Minimality Principle, Corollary 2.2, in 

m)- 

It turns out that the sequence (Viv) is uniformly bounded on compact 
sets of B. Indeed, fix i?i < 1. There exists a solution U blowing up on the 
boundary of the ball of radius 1 and satisfying U > v_m Br^ , see Remark 
2.9 in [6]. By minimality, v_ < Vn < U in Br^ , whence (Vn) is uniformly 
bounded on for any given R2 < Ri- 

We have just proved that each Vn is radial and that the sequence (Vjv) 
is nondecreasing and bounded on compact subsets of B. By standard 
elliptic regularity, it follows that (Vn) converges to a radial solution V of 
(HJ, such that V >vm B. □ 

Lemma 2.4 Assume ([3]) and (|2ip hold. Letu be a solution of ([l}. Then, 
there exist two radial functions U, V solving such that 

U <u<V in B. 



Proof. Let —M denote the minimum value of u and let v denote the 
subsolution given by Lemma 12.21 In particular, < it. By Lemma 12.31 
there exists a solution f/ > u of ([T]) and we may asssume that U is the 
minimal solution relative to v i.e. given any other solution u > v oi Q, 
U < u. In particular, U < u. It remains to construct a radial solution V 
of ([1]) such that u <V. To do so, we fix _R < 1. By Lemma 12.31 letting 
v=u\g , there exists a radial solution v = Vr of 



(24) 



Av — f{v) in Br, 

V = -I-CX3 on OBr, 
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such that Vr > « in Br. Since Vr is constructed as the monotone limit of 
minimal solutions Vn (see the proof of the previous lemma), one can easily 
check that the mapping R i-^ Vr is nonincreasing (hence automatically 
bounded on compact sets of B). Hence, as 7? — >■ 1, Vr converges to a 
solution 1/ of (HJ, which is radial and satisfies V > u in _B, as desired. □ 



3 Asymptotics of radial solutions 

Our next result establishes that the asymptotic expansion of a radial blow- 
up solution is unique. More precisely, consider the one-dimensional prob- 
lem 

(25) ^ = /(0), r<l, (^(r)^+«D as 1". 

All solutions are given implicitly by 

/ — , — 1 — r, where F' = f. 

./2F{s) 

We recall the following fact, first observed by C. Bandle and M. Marcus 
in [2] : 

Remark 3.1 Let 4> and (f)c denote two solutions of (|25|) corresponding to 
the antiderivatives F and F + c, respectively. Then <j}{r) — (t>c{r) — >■ as 
r ^ 1". 

We improve this result in the following way. 

Theorem 3.2 Let N > 1 and let ui,U2 denote two strictly increasing 
functions solving 



(26) 



d^u N - 1 du ^, , 
dr'' r dr 
lim u{r) = -I-CX3. 



Then, 

'■+°° dt 



ui{r) - U2{r)\ < C dt. 



U2(r) -P(i) 

In addition, the quantity \F{u-i) — F{u2)\ is hounded. 

Remark 3.3 Clearly, Theorem \l.l\ follows as a direct consequence of Re- 
mark WnX Lemma \2.4\ and Theorem \3.^ 

Proof. We want to think of the second term on the left-hand side of 
equation (|26[l as a lower order perturbation as r — !> 1. So, we integrate 
(|26[l in the same way we would solve (|25[). namely we let v = du/dr and 
multiply the equation by v. We get 



d fv^\ N -1 2 d 

— — H -v = — iF(u)) 

dr \ 2 J r dr^ ^ " 
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We define the resulting error term by 

(27) g:=-f.+F{u), 

which, seen as function of r, satisfies the differential equation 

(28) -i^ = v^. 

dr r 

Since m is a strictly increasing function, the change of independent variable 
u = u{r) is valid. Thinking of g as a function of the variable it, we have 

d9 = dgdr ^ Ids 
au dr au v dr 

(29) ^ = 

du r 

Since (|26p holds for r close to 1, the above equation holds for « in a 
neighborhood of +oo. Solving (|27|) for v, we finally obtain 



(30) 



/ = « = \/2 J'H-ff , 

du r r 

We start by calculating the leading asymptotic behaviour of g at +cxd : 
Lemma 3.4 

lim |M=o. 

In addition, 

(31) lim , ,1-,, , ~ 1, where G is any antiderivative of V2F. 
^ > „_>+oo (A/^ - l)G(u) ^ ■' 

Proof. First, we claim that 

(32) lim #4=0. 

u— >+oo Fyu) 



Indeed, fix e > and recalling that ([Sjl holds, choose Af > so large that 

r + oo 
Jm 

that 



f , < e. By the definition of G. there exists a constant Ca/ such 



G(ii) = Cm + / y/2F{t)dt. 

J M 

Since -F is nondecreasing it follows that 

G(m) < Cm + 2F{u) / < Cm + 2eF{u). 

JM ^j2F{t) 

Dividing by F{u) and letting u — >■ +cxj, (|32|l follows. Next, we claim that 
(33) lim 1^=0. 

u-> + oo t (U) 
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Note that by l|29p . g{u) is increasing, thus it is bounded below by a con- 
stant c as u — >■ -I-CX3. Hence, by (|30p . 

TV- 1 



^ < '-y^V^Fiu) ~c)< 2{N - 1)\/2(FH - c), 

where the last inequality holds if r > 1/2 i.e. if u is sufficiently large. 
Integrating on a given interval {uo,u), we obtain 

c < .g(^i) < 9{uo) + 2{N - 1) / \/2(F(t) - c) dt. 
Using (l32l) and the fact that limt^+oo s/^^W ^ _ i deduce (l33l) . 



Now that (|33p has been established, we return to (|30p and infer that 
given e > 0, we have for sufficiently large u, 

^ > V2(l - e)F{u) >{N- 1) V2(l - e)F{u) 

au r 

and 

g < +s)Fiu) < ^^2{1 + e)F{u). 

Integrating the above, we finally obtain for large u, 

(l_e)(7V„i) f ^2F{t) dt < g{u)-g{uo) < (l+ef^^{N-l) f y'2F(t) dt 

and H3ip follows. The fact that g{u) — +cxd as u — > +oo follows automat- 
ically. □ 

Next, we prove that given two solutions ui,U2, the corresponding error 
terms g\ , g2 given by ((27| differ by a bounded quantity. 

Lemma 3.5 Let iti and U2 be two solutions of (|26p . Introduce Vi — 
and ^ 

g^ = -Y+^M^ /or i = 1,2. 

Then, gi — g2 is bounded. 

Proof. We begin by rewriting the system (|30p as a nonlinear integral 
equation with unknown g. To do so, solve the first line of H30p for r: 



Differentiate with respect to u: 



du 



2 



Equate the above equation with the second line of pop . to obtain the 
following second order differential equation: 

rf'ff , 1 1 (dg\^ f-t dg 



(34) nil -I : : — —] — 21;- -^ = o 

^ ' du'2 2{N - 1) F - g \du J 2{F-g)du 
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Define 

So, 1/{F — g) — 1/F + q and p4|) can be rewritten as 



+ ■ 



+ q 



du'^ 2(iV-l) \F ' \duj 2 \F 

In other words, k = dg/du solves the logistic equation 

N 



fc - / = 



The general solution of such an equation is well-known and is given by 
2(iV-l) 



k 



N 



where uo,C are arbitrary constants. Since all integrands are positive for 
large u, we may take ito — +oo in the above formula and obtain 



(36) 



k = - 



2(N- 1) 



iV 



i:^{^+q)e-i-r^"^''f'^dt + C 



2{N - 



N 
(iV-1) 



^-Vf- 



e 2 ■ 



N 



Next, we identify the leading asymptotics of the quantity /^°° ifdt. To 
do so, simply recall the definition of q given by (|35[) . as well as the leading 
asymptotics of g given by Lemma [3.41 

(37) I qfdt^l ^^^^ 

where G' = V2F. Integrating by parts, we discover that 



(38) / ^Cl^dt 



G /■+°° V2F , 

h / dt ■ 

F L. F 



Using this in p6|) . we deduce that 

iV- 1 1 



+ 2 



dt = 0(1) 



N C 



2F. 



In addition, k = dg/du = {N - l)/r ^2{F - g) {N - 1)V2F. So, we 
must have C — —1/N and so 
(39) 

(iV-1)^ 



du 



N 



'2F 
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Take now two solutions ui , U2 of (|26|) and let gi , g2 denote the associated 
error terms. By (|39|l . we have 



du 



du 



N - 



N 



'2Fx 



r + oa / 1 



^/2F ^ 2 V y e iV 



r+oo /' 1 



where q = q, satisfies (|35|l for g = Qi- Reducing to the same denominator 
and using (|37l), ((Ml), it follows that 



dgi dg2 
du du 



N- 



7V3 



2F 



■ /+~ nfdt 



e 



1 +fV2Fle-S 
/2F 2 



7^ + ^ 
/2F 2 



iV 

1 " 

iV 



/2F 1 e-i^^*°°''^^'''dt 



To simplify the above expression, we write Ci = e 2 /u*^ 9;/^* -yy-g obtain 



rfgi _ rf.g2 

du du 



ei 



62 









/'+00 


62 










TV- 1 



2F ( 61 
1 



92 



/2F 



+ ■^V2F 62dt- — 



A/' 



+ ^V2F\eidt^ 

2F 2 y 

TV- 1 



TV4 



/2F(6i -62) 



ATS 



-L + f ^/2F)62dt 

/2F 2 



- 62 



1 

" TV 
'2Fx 
1 

v/2F 



91, 



2F eidt 



TV - 1 ; — , , TV - 1 ; — 

^(61 - 62) tt:^V2Fx 



(61 — 62) 



ATS 

= + ^V2F] eidt 
2F 2 



^ + SlV2F]e,dt 
/2F 2 / 



-L + ^V2F ) 62df 

/2F 2 



Cancelling lower order terms in the above expression and noting that 
62 ~ 1, we obtain 



r-\-oc 



^ + ^V~2F]e^dt 
/2F 2 



4-00 



^ + f ^/2F ) 62dt 
/2F 2 
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The right-hand side in the above expression can be rewritten as 



iV3 



7V4 



1 1 f+°° — 

^{ei—e2)dt+- / V2_F (giei — 9262) 
2F 2 Ju 

V2F {qiei - (7262) dt 
2F(gi(ei-e2)+e2(gi-g2))dt 



^ ^ /2F(ei -62) - ^^ttt!-\/2F 



^-^/2i?(e.-e2)-^^ 



Since ~ 1, we have, using the mean value formula. 



(40) ei - 62 

In addition, by (|35|l . 



(<?! - q2)fdt. 



qi ^ and 51-172 



F2 ■ 



So, 



ei - 62 ~ - - 



and it follows that 



du du 
AT- 1 



2Ar3 
Hence, 



1 Af- 1 

2 iV* 

+ 00 

2F / ^/2F 



^2F 



f 

71 -g2)-^dt- 



7 



F2 



dt 



dffi dg2 


<Cv^(/_ 


du du 





G 

p3/2 

< CV2F 



\gi -92\^ dt+ 



\gi-g2\-^dsdt + 



m - 92\ 



dt 



1.91 - .92I dt- 



(2F)3/2 

\g1-g2 



(2F)3/^ 



where G' = y/2F. We want to estimate further each of the two terms 
on the right-hand side of the above inequality. Since gi = 0{G), one can 
easily check that all integrals are convergent. In particular, we may always 
find [/ > It so large that 



f f'^ f 

\gi-g2\-^dt< J \gi-g2\ — dt, 



\gi - gal 
(2F)3/2 



dt < 



|gi - gal 

(2F)3/2 



dt. 
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It follows that 

<CI sup Isi-jjII |-i= + v^ * 



te[u,u] 



/2F A (2F)3/ 



< sup jpi -32I 

Integrating the above expression between a given constant uo and u, we 
obtain 

" dt 



\gi - 92\(u) < \gi - g2\{uo) + C \ sup I51 - 32 

Choose now uo so large that C f^°° < 1/2. It follows that 
sup i.gi - 32I < 2lgi - 32|(tio) = Co. 

teluo,U] 

This begin true for U arbitrarily large, we finally deduce that gi — 32 is 
bounded, as desired. □ 



= 1. 



Completion of the proof of Theorem 13.21 

Let u\,U2 denote two solutions of \2Q\ . By pO[) . each ut, i = 1,2, 
solves 

dui I dr 
^2{F{u,) - g,) 

Integrating, we obtain 
^+00 ^ 

dt = l- r ^ I dt. 



1 ^2{F{t)^gi) Ju, y/2{F{t) - 32 

Without loss of generality, for a given r we may assume U2{r) > ui{r). 
Split the left-hand side integral : f^°° = f"^ + f^°°. It follows that 



r^2 H 
Jm ^2{F{t)-gi) ^ 



^2{F(t)-g2) ^2{F{t)^g,) 



dt 



, ^2{F{t)-gi)^2{F(t)^g2) 



3^-3' dt 



^2{F{t) - gi)^/2{F{t) - g2) (v/2(i^(f) - gi) + ^/2{F(t) - 32)) 

Recall that by Lemma [3.4l gi — o{F) as t — > +00. Recall also that 32 — (?i 
is bounded. So, for sufficiently large values of U2, we deduce 

(41) / ^ ^ < C ' 



1 v/2F(I) - Ju, F{tf/ 
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Since F is increasing, it follows that 



< ^^1=^ < C I dt<C 



1 , , „ /■+°° dt 



C f+°° dt 



< 



Hence, 

r+°° (It 

as stated in Theorem 13.21 It remains to prove Witliout loss of gener- 
ality, we assume ui{r) < U2{r) so 

dt /■+°° dt /■+°° dt 



/■+°° dt 



3 ^F(t-(u2-wi)) 

+°° VF(t-(u2-m)) 

- (1*2 - til)) 

F(t)-F(t-(i.2-m)) 

2 ^jF{t)F{t - (it2 - ui)) (/FW + - ("2 - u^))) 

" dt 



Recalling (gl]), (O follows. □ 



4 Gradient estimates 

Proof of Theorem 11.21 Let w = ui — U2 denote the difference of 
two solutions. Without loss of generality, we may assume that it2 is the 
minimal solution of ([1]), so that ui > U2 and U2 is radial. 
Step 1 : estimate of tangential derivatives 

We begin by proving that any tangential derivative of w is bounded. 
Since the problem is invariant under rotation and since U2 is radial, we 
need only show that (r, 0, . . . , 0) remains bounded as r — >■ 1~ . Given 
X = {x\,X2,x') £ B and 6 > Q small, we denote by xe = (a;icos^ — 
X2 sin 6, xi sin 9 + X2 cos 6, x') the image of x under the rotation of angle 6 
above the xi-axis in the (xi, X2) plane. By the rotation invariance of the 
Laplace operator, the function ue defined for a:: £ B by 110(2:) = ui{xe), 
solves U]). Using (Q and assumption we deduce that wg = ui — ug 
solves 



Awg + Kwg = fiui) — f(ug) in B, 

wg = Q on dB. 
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By the Maximum Principle on small domains, there exists Ro € (0, 1) such 
that the operator L = A + K is coercive on i? \ Br^ . As a consequence, 
we claim that there exists a constant C > such that for all x £ B\ Bjig , 

(43) \'wg{x)\ < C sup \wg\ . 

Let indeed C > denote the solution of 

in B \ Br„ 
on dBR„ 
on dB. 

We shall prove that :— wg — ±supg^^ l^^elC respectively non- 
positive and nonnegative, which implies that (|43p holds for the constant 
C = IICIIoo. We work with z'^ and assume by contradiction that the open 
set u} — {x £ B\BRg : z'^{x) > 0} is non-empty. Restricting the analysis 
to a connected component, we have 

Az+ + Kz+ = f{ui) - f{ug) > in w 

z"^ < on duj. 

By the Maximum Principle, we conclude that z'^ < in oj, a contradiction. 
We have thus proved (|43|l . Since iti G C^{Br„), we deduce that for some 
constant C > and all x £ B \ Br^ , 

\wg{x)\ < C9. 

Applying the above inequality at the point x = (r, 0, . . . , 0), r £ {Ro, 1) 
and letting S — > 0, we finally deduce that 



<C for aU r e (i?o,l) 



as desired. 

Step 2 : estimate of the radial derivative 

It remains to control dw/dr. Fix R G (0, 1). Let GR{x,y) denote Green's 
function in the ball of radius R. Then, for x £ Br, 

w{x) = [ ^^{x,-)wda +[ GR{x,-){f{ui) - f{u2)) dy 

(44) ieSR OUy Jg^ 

^■.Wi{x) +W2{x). 

We want to let i? — >■ 1 in the above identity. To do so, we first observe 
that wi is harmonic. By the Maximum Principle, |wi| < ||w|jicxj{ga^). 
By estimate (|4]), we conclude that uii — >■ as ii — ^ 1. To estimate W2, we 
need the following crucial estimate : 

Lemma 4.1 Assume (|6}. Then, 

sup / l/(ui) - /(U2)l (r,6l) dr < +00. 
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We shall also need the following elementary estimates. 

Lemma 4.2 There exists a constant C > such that for all 1/2 < r, R < 
1 and all x,y £ Br, 



(45) 



(46) 



Gi(x,-) da < 1 



dGi, , 



da <C 



We postpone the proofs of the above two lemmas and return to (|44 
Using polar coordinates, 



W2(x)= GR{x,-){f{ui) ~ f{u2}) dy 

J Br 

L {Lb ^^(^■'■)(/("i)-/(^^2))rfa^ dr 



By Lemmas 14.11 and 14.21 we may easily pass to the limit in the above 
expression as _R — >■ 1, so 

w{x)= f Giix,-)if{m)~ f{u2))dy 

J B 



Using again Lemmas 14.11 and 14.21 we also have that w is differentiable in 
the r — \x\ variable and 



Using polar coordinates again and Lemmas 14. II and 14.21 we finally obtain 

dw 



< C+ 



dr 



+ sup 

r6{l/2,l) \JdBr 



dGi 



da] sup / \{f{ui)-fiu2)\{r,9)dr 

< C. 



It only remains to prove Lemmas 14.11 and 14.21 

Proof of Lemma 14. 11 We first deal with the case where ui,U2 are radial 
and iti > U2. By assumption (|6]), we have 

\f{ui)- f{u2)\ dr<J^ (/>i)-/(u2)) dr + A-||«i -U2||£oo(B). 

Using Q, we see that ui — U2 is bounded and so it remains to estimate 
/(wi) — f{u2). By pO|l. each ut, i — 1,2, solves 
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We also know by Lemma [3.41 that gi = o{F{ui)). So, 

,. dui/dr 

lim — , — 1. 

Using this fact, as well as Lemma [3. 51 and ((5|, we obtain for R £ (1/2, 1), 



/' 

Jo 



R _ rR 

(/(wi)-/(M2))dr< / {f{u^)- f{U2)) dr + K\\u^-U2\\L^(^B) 

Jo 

< C Hfiu,) - fiu,))^^^ dr + C 
Jo \/2F{u\) 

<C ri nu,)^^^ll^ - fiu,)J^^] dr+ 
- Jo \ W^F{^) '^/2F{^J 

+ C rfiu,)(^^0^-^^^^] dr + C 

< c (\/FK) - /fm) {R) + C+ 

Jo V V^Fi^ J 

+ C j(u2)— , = dr 

Jo y/F{u2)F{m) 

VFiMR)) 

+ C r /(^^) giF(u2)~g2F{m) 

Jo y/F{u2)F{ui) y^F{u2)F{ui) 

<C + C r /(^^) {gi-g2)F(u2)+g2iF{u2)-F(m)) 
Jo ^F{u2)F{ui) ^F{u2)F{ui) 

<C + C||ffi-52||^oo(fl) r iipldr+ 

'Jo F{U2) 

+ C\\F{u,)~F{u2)\\l^(b) I -J^dr 

Jo ^(U2) 

<C + C r l^^J0^dr + C 

Jl/2 F{U2) ^2F{U2) 

<C + C (^F-^''^(l/2) - F"^/^(i?)j < C. 

This proves the lemma for radial solutions. To obtain the estimate in the 
general case, we may always assume that U2 is the minimal solution of 
([TJ, so that U2 < ui and U2 is radial. By Lemma [2.4l up to replacing / by 
/ given by Remark 12.11 there exists another radial solution V such that 
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V > ui > U2. Using assumption (|6]), we have 

- /(u2)i dr< j (/(?ti) - /(U2)) dr + A'llui -U2||L~{i3) 

' [KV) ~ f{u2)) dr + K\\m - U2|Uoo(s) 

< ifiV) - f{u2)) dr + 2K\\m - U2h^(B) 
Jo 

By Q , 111 — U2 is bounded and the result follows from the radial case. □ 

Proof of Lemma l4.2l (|45p is standard : write the representation formula 
(|44|) both in Br and in Bi , change variables in the Bi integral and identify 
the kernels. Next, we prove that given any r £ (0, 1), Jg^ Giix, ■)da < 1. 
It suffices to show that for any (j> G Cc(0, 1), 

(47) ly^^'^^laB ^^(^'■)'^'^) ^'"^^ ll-^lli^o.i)- 

By definition of Green's function, the left-hand side of the above inequality 
is the function v solving 

-A« = <j> in B, 

V = Q on dB. 

The above equation can also be integrated directly : 

i/(r) = r^-^ /' mt""'^ dt, 
Jo 

whence < ||(/>||i,i(o,i) s-nd \v\ < ||i/>||i:,i (o,i) i-e. (f47l) holds. This proves 
that Jgg Gi{x,-)da<l. 

We turn to the second estimate in (|46[) . Recall that the Green's func- 
tion in the unit ball is expressed for x,y £ B, x ^ y,hj 
(48) 

Gi{x,y) = r (^{B? +r'^ -2Rrcosipf^'^^-V (^(l + - 2Rrcosipf''^^ , 

where R= \x\, r — \y\, ip is the angle formed by the vectors x and y and F 
is the fundamental solution of the Laplace operator. Differentiating with 
respect to R, we obtain for some Cjv > 0, 

(49) CMf^^{x,y) = 

R — r cos ip Rr^ — r cos ip 



+ r2 - 2Rr cos ip^/"^ (1 + R'^r^ - 2Rr cos ip)^/'^ 
_R — r 4- r(l — cos (jp) _Rr^ — r + r(l — cos oj) , „ 



{{R - r-)2 + 2Rr{l - cosip))^^^ ((1 - Rr^ + 2Rr{l - cos ip))^^^ 
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We estimate A and leave the reader perform similar calculations for B. 
Clearly, given e > 0, the expression (|49p remains uniformly bounded in 
the range 1/2 < R,r < 1, e < ip < 2-k — e. Hence, 



\A\ da<C, + C / \A\da. 

For y € dBr n [0 < <^ < e], let z = z{y) denote the intersection of the 
line (Oy) and the hyperplane P passing through x and tangent to the 
hypersphere dBn. Then, there exists constants ci,C2 > such that for 
all y G dBr n[0 <(p <e], 

ci (1 — cos (j)) < \z — <C2(1 — cos 0). 

Hence, letting B^~^{x, p) C P denote the A''— 1-dimensional ball of radius 
p > centered at s, we obtain 



^ I , ^ ^ I , f \R — r\ + Cr \z — x\^ 
A\ do- < C 1 + / ^ ^ dz 

9Br \ JB'^--l-(x.RsinE) [\R-r\^ +c\z- x\ ) 

<c{i+i '^'T' + '^i!!' dz 

Ibn-i(o,r,) {\R-r\' +c\z\') ' 
<C(1+ / \R^r\+C\R-r\^\z\^ 



<C\l+l ^-—^dz + \R~r\ I U?-"" dz 

< C. 



Working similarly with the B term in (|49[). we finally obtain the desired 
estimate (Hell. □ 



5 Asymptotics to all orders 



This section is devoted to the proof of Theorem II. 9 1 Our first task consists 
in applying the Fixed Point Theorem to the functional N defined for 
V € B{vo, p), u € I hy 



(50) W{v)]{u) = y^2 (^F{u)-{N-1) £ld?j, 

where r is given by (|13|l . Let us check first that M {B{vo, p)) C B{vo,p). 
Take v G B{vo,p). Then, 

(51) 1 > r > 1 / — dt=l / ^= dt. 

By ((3)1, it follows that for p < 1/4 and Uo sufficiently large, 1 > r > 1/2. 
Hence, 

/ - dt <C V2F dt = o{F(u)), 
Juo ^ Jua 
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where we used Lemma 13.41 So for Uo large and u> Uo, 
N -1 



- dt 



Uo 



We deduce that 

N(v) - Vo 



(52) 



N ■ 



F{n) 



1 r V 
) Juo 



dt < 



N ■ 



F{u) 



-dt 



< P- 



Next, we prove that N is contractive. Given V\,V2 £ B{vo,p), let ri = 
r(u,vi),r2 = r{u,V2) (where r is given by (fT3|l ). Then, by estimate (fSTj) . 
1/2 < ri,r2 < 1 and 



AA(«i)-AA(«2) 



DO 



1 - 



N_-l_ 



Ua 



^dt- 
ri 



N ~ 



F{n) 



1 r «2 

) Juo ^2 



dt 



< C 



N - 



F{u 
C 







V2 


Iuq 


ri 


r2 



F{u) 

C 
F{u) 

C 
F(n) 



Cp 
F(u) 

Cp 
F{u) 



Uo 



Uo 



dt 

\vi — V2\ dt + 

u 

V2F dt + 
2F dt+ 



r 


1 


1 


dt^ 


/ vo 






JUo 


ri 


r2 





'2F\ri -r-2| dt 



Uo 



Uo 



Uo 



'2F 



I2F dt + 



1 

V2 



Uo 



I2F 



Uo 



ds 

+ 00 



ds 



dt 



'2F dt. 



Uo 



Using Lemma [3.4l we conclude that M is contractive in B{vo, p) if Uo was 
chosen large enough in the first place. We may thus apply the fixed point 
theorem. 

So, it only remains to prove (|f 6p . We first observe that the sequence 
(wfc) defined by (fT4|l is asymptotic i.e. Wfc+i(M) = Vh{u){l + o(l)), as 
u — >• +0O. Since Hfc+i = Af{vk), it suffices to prove that Af{vo) — vo = o{vo) 
and iterate. By (|52[) . 



A/"(tJo) — VQ 



< 



C 
F{u) 



2F dt 



Uo 



and the claim follows by Lemma [3.41 So, the sequence (vk) is asymptotic 
and so must be the sequence (it^) defined by (|f 5|l . We are now in a position 
to prove (|f 6|l . By Theorem If .11 we may restrict to the case where u is 
radially symmetric. Let v = du/dr. By (|26|) . v solves 



cfo ^ N - 
dr r 



-V = /(w) 



costin-dupaigne-jdif f eq-janlO . tex 



March 19, 2010 



25 



Use the change of variable u = u{r) to get 

dv N-1 

V- 1 V = f{u). 

au r 

Integrating, it follows that for some constant C 

(53) ^ = F{u) + C- -V dt. 

2 Juo r 

Up to replacing F{u) by F{u) — F{u) + C (which is harmless from the 
point of view of asymptotics) , we may assume C = 0. So it suffices to 
prove that v € B{vq, p) to conclude that v coincides with the unique fixed 
point of A/", whence (|16|) will follow. By (|53|) (with C = 0), v < vo and so 

0<Vo-v< y/2F{u) - ^2 (^F{u) ~ /" ^ ~ ^ 
f" y/2Fdt 

~ V^fT 

By Lemma [3.41 it follows that 



Vo — V 

< < Vo ~ V < p 

vo 

and V € B{vo, p) as desired. □ 



6 Universal blow-up rate 

In this section, we prove Theorem II. 101 that is we characterize nonlinear- 
ities for which the blow-up rate is universal. 

Proof of Theorem [TTTni . 

Step 1. We begin by establishing the theorem when Q — B is the 
unit ball. In light of Theorem 11.11 it suffices to prove p9[) for one given 
solution u of IT]), which we may therefore assume to be radial. By (|30p . 
we have after integration that 



(54) / , dt 

Ju V^{F{t) ~ g) 

By definition of uo, we also have 



(55) / , dt = 1 - r. 

Observe that u > uq, split the integral in (|55p as /^^°° ~ J^^ + J^°° a-nd 
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equate (|54p and (|55|) . It follows that 



dt= — . dt 



y^mt) Ju W2iF{t)-g) y^2F{t) 



f 



y2F(7))\/2(F(t) - ») (\/2F(I) + v/2(F(f) - g)) 



dt 



Recall that by Lemma r3.4l g = o{F) as t — >■ +oo and g{u) ~ (A*'— l)G(u) 
(A'' — 1) \PTF dt. So, for sufficiently large values of u, we deduce 

56 / , dt<C ,,,,, dt. 

Since F is nondecreasing, it follows that 

Uq 



< , < / , dt<c ,,,,, dt. 



Hence, 



r + oa r* 

< tt - ito < Cv/2F(tt) / f- 
J u (2 



2F ds 

dt. 



(2F(t))3/2 



and ((191) follows from (fT8|) . 

Step 2. Next, we prove that (|19p holds for general domains fi. To this 
end, we combine a standard approximation argument by inner and outer 
spheres (see e.g. [9]) and the comparison technique of [6]. Let u denote a 
solution of (|17|) and take a point xo G 90. Let B G Q denote a ball which 
is tangent to d^l at xq. Shrink B somewhat by letting B^ = {1 — e)B, 
e > 0. Observe that u £ C{B^) is a subsolution of 



(57) 



Af/ = /([/) in Be 

f/ = +00 on 



By Lemma [2.3l there exists a solution Ve of (jSTp . such that > u m B^. 
Furthermore, 14 can be chosen to be the minimal solution of (|57p such that 
Ve > u in Be. In particular, Vs is radial and £ — >■ 14 is nondecreasing. In 
addition, e — >■ 14 is uniformly bounded on compact subsets of B (working 
as in the proof of Lemma l2.3p . so 14 converges as £ — >■ 0, to a solution V 
of H]) such that V > u in B. By Step 1, 

i}i^^V(x)-uo{l-dB{x)) = 0, 

where ds denotes the distance to dB. Since V > u and since the above 
discussion is valid for any point xo £ d^l, we finally obtain 

(58) limsup[u(x-)-tto(l-d(a:))] <0, 

x^an 
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where d(x) is the distance to dfl. Choose now an exterior ball B C \ 
which is tangent to dfl at xq. For e > small and _R > large, the annulus 
— RB \ (1 — e)B contains Q. Let Ue denote a large solution on A^, 
which we may assume to be minimal, radial and bounded above on Q 
by u. Again L^^ — >■ C/ as e — >■ where f/ is a radial large solution in 
A = RB \ B D Repeating the analysis of Step 1. (which was purely 
local) for the case of a radial solution defined on an annulus rather than 
a ball, we easily deduce that 

\im U{x) - uo(l - dsix)) = 0. 

Since u > U and since the above discussion is valid for any point xq £ d^l, 
we obtain 

(59) liminf [u(x-) - uo{l - d{x))] > 0. 

So, by (|59p and (|58|) . we have that (|19|) holds in any smoothly bounded 
domain SI. 

Step 3. It only remains to prove that H19p fails when (|20p holds. We 
use Theorem ll.9l to compute the second term in the asymptotic expansion 
of a solution. By (fT4)) . 



whence 



vi{u) = ^l2{ F{u) - (N - 1) V2F dt{l + o{l)) 



1- ' ^ + ,„-„JI^,i+„„, 



1 r V2F dt 

+ (^-l)7^^73Ti72(l + «(!))■ 



Integrating (|15p for fc = 1, it follows that for r close enough to 1, 
r + oo 1, r + oo r* ./o p Jo 



Recall (|55|l . split the integral in (|55|l as = Iu°° ^''^'^ equate 

(|60)l and ([551) to get 



Since F is nondecreasing, we deduce that 
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Note also that 



(62) 



"1 J*V2Fds 
(2F)3/2 



< / — dt< 
2F 



Assume by contradiction that hm, >i- (iti— ito)(r) — 0. Then, (|62|l implies 

that 



(2F)3/2 



\y2F(V 



Using this information in (|6ip . we obtain that 



Ml — Uo 



/■+00 f* ^/2F ds 



But (|20p would then lead us to a contradiction with the assumption 
lim^^i- {ui — uo){r) = 0. So we must have 

liminf (ui — uq) > 0. 

r— 



and so p9[) fails. 



□ 



7 The first three singular terms 

In the previous section, we characterized nonlinearities for which only one 
term in the expansion is singular. In the present section, we calculate 
implicitly the next two terms in the expansion. We have not tried to 
characterize those / for which all remaining terms are nonsingular, but 
this can certainly be achieved. We leave the tenacious reader try her/his 
hand at this computational problem. 

We begin by calculating the leading asymptotics of vi, V2. By (|14p . 
we have 

2 ru 

^^F-{N-1) J V2F{l + o{l))dt. 

So, 

^=V^(l-(iV-l)£^^(l + o(l))' 



In other words, 



^2F 



To calculate V2, we introduce some notation. Given a positive measurable 
function v, set 

Fv = / V dt, Qv = , iiv — I — , 

./ V ./„ V 
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and 

Tv = {N - l)PQv + P{vRv). 
vi is then expressed by 

vi=vo-{N-l){l + o{l))Qvo, 

while V2 is given by 

^i=F-{N-i) r — dt= 

= F - (TV - 1) J {vo -{N- l)Qvo + o{Q{vo)){l - Rvo + o{Rvo)) dt - 
= F-(N- l)Pvo + (iV - l)Tvo{l + o(l)). 

So, 

V2 = (2F - 2{N - l)Pvo + 2{N - l)Tvo{l + o(l)))'/'' = 

= ^0 (^l-{N-l)^ + iN- 1)^(1 + o(l))) = 
, N-l PvQ , N-l TvQ 3 .^2fPvo\\ 



+ 



o{Tvo/F + iPvo/Ff)y 



And so, 



i;2fo\ 2 F 2 F 8 \ F J 

+ o{Tvo/F+{Pvo/Ff)^ = 

^ + - 1)^ + - 1) (-^ + - 1)^1 (1 + oil)) 



1 /^r .J'-^dt 



2F ' ' (2F)3/2 
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Integrating once more, we finally obtain 



x{N-l) (- i{N-l) 

Ju^(r) \ J \ 




This proves Proposition 1 1.141 

8 An example: f(u) = u^, p > 1 

Finding the 71-th term in the expansion for abitrary 71 G N is out of reach 
for general /, simply because of the algorithmic complexity of calculations. 
However, when additional information on / is available, one can guess the 
general form of the expansion and then try to establish it. This is precisely 
what we do in this section, with the nonlinearity f{u) = u'' , p > 1. 

For notational convenience, we shall work with F{u) = ^u^'' , where 
2q — 1 = p, which simply amounts to working with a constant multiple of 
the original solution. 

Recall (|14p and (|15p . We want to prove inductively that there exists 
numbers ak, bt depending on k,p, N only such that 

n 

(63) v„ = u''J2 bkU'"^"-'^ + o(^.«-"(«-^'), 

n 

(64) it„ = ^afed'-' + £'„(d"^ 



where £;„(d"~+"+^) ~ e„d —+"+1 for some e„ G K, as d 0+. We 
have vo = y/2F = u'' . Solving for uo in (fT5)) yields uq — c d~'^ . So, 
and l|63|l hold for n = 0. Suppose now the result is true for a given n € N. 
In the computations below, the letter denotes a number depending on 
k,p,N only, which value may change from line to line. By (|63p . we have 



\fe=0 
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So, 

(65) r°°£f£=ti-«(^^c,t-^(«-^A +o(f-("+i)(''"i)) 

(66) = Xi^**"'^'"" + 



It follows that 



, _ r + oo ds_ 
■It v„ 



(n + l)(<7-l)N 



Whence, 



And so, 



v„{t) 



1 _ r+oo ds_ 



= j2F-(iV-l) / -—^^dt = 



fc = 

n + 1 

l + ^c,K-'=(«-^)+o(u-("+i)(«-^') 

fc=l 

n + 1 



fe=0 

This proves (|63p . Integrating (|15p . we obtain 

(67) / - = d. 

Now, =i;„+c„+im''-("+i'(''-i' (1 + 0(1)). So, 

^ = ^ + c„+i«-^-("+^»'-^) (1 + o(l)). 

It follows that 

/' + °° du /' + °° dtt _(„ + 2)(9-l) , f^ss 

d= — +c„+ni„|i '(l + o(l)). 

In addition, u„ ~ vq, so u„ ~ uo, and so ~ d. Using this in the 

above equation, we get 

d+cn+id"+\i+o{i))^ ^, 
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Recalling that Vn is defined by (|67|) and satisfies (|64|l by induction hy- 
pothesis, we conclude that 

1 " 

Vr, + 1 ={d + c„ + id" + 2(l + o(l))"~ ^ afe (d + c„+id" + 2(l + 0(1))' 
Expanding again the above expression, we finally obtain 

n + l 

V„ + l = d-T^ "■"'^'^ + Sn + l(d"^+" + '), 



which proves (|64p . Proposition 1 1 . 1 5l follows . 
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